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In  this  report  we  derive  the  solution  of  the  3f  -minimal  weighted 
sensitivity  problem  with  1  pole/1  zero  sensitivity  weighting  functions 
for  stable  minimum  phase  rational  plants  with  an  input  delay.  We  calcu¬ 
late  the  explicit  feedback  compensator  and  analyze  its  stability. 
Although  we  have  a  technique  for  approximating  these  compensators  with 
proper  compensators,  the  full  description  of  this,  along  with  an  analy¬ 
sis  of  stability  robustness  will  appear  in  a  future  report.  Another 
report  will  also  cover  generalizations  of  this  work  to  plants  which  are 
unstable,  to  plants  which  have  right  half  plane  zeros,  and  to  the  case 
of  a  more  general  sensitivity  weighting  fiinction. 

Preface. 

The  results  presented  here  expand  upon  the  work  Initiated  in  Appen¬ 
dix  B  of  [Flamm  1985a]. 


Purpose  of  Research. 


The  goal  of  this  work  is  to  obtain  euid  analyze  explicit  compensa¬ 
tors  for  delay  systems  which  achieve  or  approximate  a  closed  loop  sensi- 

99 

tlvity  function  minimal  in  the  H  -norm.  A  major  premise  of  this  inves¬ 
tigation  is  that  all  real  systems  contain  delays,  so  that  it  is 
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only  by  examining  the  optimum  for  such  systems  that  an  understanding  of 
achievable  performance  can  be  obtained. 

00 

We  shall  also  examine  the  robustness  of  the  stability  of  3f  -optimal 
sensitivity  compensators  with  respect  to  changes  in  the  delay  present  in 

CD 

the  plant.  We  show  below  with  a  simple  example  that  1  -minimal  sensiti¬ 
vity  feedback  systems  are  not  necessarily  robustly  stable  with  respect 

to  the  addition  of  a  small  delay  in  the  feedback  loop.  This  means  that 

00 

there  is  something  wrong  with  the  W  sensitivity  problem  ais  commonly 
treated.  As  a  consequence  of  this  lack  of  robustness,  the  resulting 
feedback  system  is  ill-posed  in  the  sense  of  [Willems  1971]  pp.  90-91. 

We  remark  here  that  the  i 1 1-posedness  for  the  feedback  system  seems 
to  be  a  result  of  the  memoryless  part  of  the  opened  feedback  loop  having 
gain  which  is  too  large.  See  [Willems  1971]  p.  100.  From  this  point  of 
view  a  strictly  proper  approximation  to  the  optimal  compensator  for  an 
3f  problem  for  a  rational  transfer  function  plemt  will  likely  result  in 
a  well-posed  system.  Nonetheless,  it  is  troublesome  to  use  the  formula¬ 
tion  of  a  problem  for  which  the  ideal  solution  is  ill-posed. 

For  plants  with  a  delay  in  the  input,  which  we  shall  consider  in 
this  report,  the  optimal  feedback  systems  will  be  well-posed,  according 
to  [Willems  1971]  Corollary  4.1.1,  precisely  because  of  the  delay  in  the 
pltuit.  (This  conclusion  depends  upon  the  nature  of  the  optimal  compen¬ 
sators  for  these  systems.  We  examine  this  later  for  the  case  treated  in 
this  report,  and  see  that,  although  these  compensators  are  unstable,  and 
in  fact  have  infinitely  many  poles  in  the  right  half  plane,  they  are  not 
anticipative. ) 

We  also  intend  to  test  the  conjecture  that  designs  for  systems  con 
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cainlng  delays  will  provide  a  means  of  obtaining  compensators  which  are 
robust  for  high  frequency  phase  uncertainty,  such  as  that  resulting  from 
right  half  plane  zeros  of  the  plant  for  which  the  designer  of  the  com¬ 
pensator  is  uncertain  of  the  location. 

We  note  that  in  the  example  used  below  to  illustrate  the 

CD 

i  1 1-posedness  of  the  W  optimal  feedback  system  for  the  purely  rational 
plant.  [Zames  and  Francis  1983]  p.  592  SVI.A,  remark  that  the  optimal 
"Q"  parameter  is  a  lead  filter.  Since  this  has  a  certain  intuitive 
appeal  (which  those  authors  suggest),  we  would  like  to  examine  how  this 
may  extend  to  the  case  of  delay  plants.  In  particular,  we  are  led  to 
ask,  in  what  manner  can  we  look  at  the  delay  in  terms  of  right  half 
plane  zeros:  for  example,  as  Infinitely  many  zeros  at  the  point  <»? 

Finally,  we  hope  to  provide  conceptual  insight  into  the  reason 
delays  limit  the  sensitivity  of  a  compensated  system. 

We  also  regard  the  delay  problems  considered  here  as  a  first  step 
towards  considering  similar  design  issues  for  other  infinite  dimensional 
plants. 

CD 

An  Example  of  II 1-POsedness  of  the  W  -minimal  Sensitivity  Feedback 
System  for  a  Rational  Plant. 


CD 

The  Importance  of  understanding  t  -minimal  sensitivity  design  for 
delay  systems  is  emphasized  by  the  fact  that  designs  which  achieve 

CD 

W  -minimal  sensitivity  for  finite  dimensional  systems  are  not  generally 
stable  when  a  small  delay  is  added. ^  We  Illustrate  this  fact  with  a 

^This  point  was  suggested  by  Gunter  Stein. 
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simple  example. 


00 

In  [Ilcunes  and  Francis  1983,  p.  593]  Che  solution  co  the  3f  -minii-nal 
sensitivity  problem  is  computed  for  the  case  of  a  stable  plant  with  two 
right  half  plane  zeros.  We  consider  this  example  here.  Let 

be  the  plant,  where  is  stable  and  minimum  phase,  and  let 

» =  •** 

be  the  sensitivity  weighting  function.  The  optimal  "Q”  parameter  is 


Q(s)  =  Tl  -  Pte-s)  -».(bt*s)(b2-t-s).p  -1 
-  Li  W(s)(c+s)-l  (b,-s)(b,-s)  • 

where  D  and  c  are  constants  determined  by  P  and  the  b^s.  D  mav  be  posi¬ 
tive  or  nesatlve.  depending  on  whether  8  >  1  or  fl  <  1.  (We  will  not 

show  it  here,  but  0  <  1  is  the  only  interesting  case.)  We  assume  D  is 

positive,  for  0  <  1.  Then  the  optimal  feedback  compensator  is 

F  =  Q(l-PQ)-^ 


Now  suppose  that  the  true  plant  is  really  P(s)  =  e”^®Po(s),  with 
e>0.  We  check  stability  for  the  closed  loop  system  with  true  plant  and 

p 

compensator  designed  for  the  nominal  plant  P  by  computing  and 

checking  the  location  of  its  poles.  A  straightforward  computation  gives 

_E_  .  p _ _ . 

«*‘D(c-s)  ♦  [W(s){c+s)-D(c-s)] 

We  now  show  that  this  is  unstable  by  showing  that  the  denominator,  call 
it  H(s).  in  the  above  expression  has  infinitely  many  right  half  plane 
zeros. 

For  this  we  evaluate  M(s)  on  the  imaginary  axis  emd  separate  it 
into  real  and  imaginary  parts.  N(lu)  =  F(u)  *  t*G(u).  Now  we  show  chat 
F(<j)  has  only  finitely  many  real  zeros,  and  conclude  by  appealing  to 
[Pontryagin  1955],  Theorems  3  and  6,  that  D(s}  has  zeros  in  the  right 
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half  plane. 

A  compucacion  gives 

1^1 

F{«)  =  u*D[cos{eu) - g— ]  +  <i>D(p-c)sin(£u)  +  PcD[cos(£(j)-1]  +  c. 

Since  D  >  0.  this  equation  cannot  be  satisfied  for  sufficiently  large  u. 
and  we  conclude  that  F(u)  has  only  finitely  many  real  zeros.  Therefore, 
D(s)  has  zeros  In  the  right  half  plane. 


Problem  Considered. 

Oft 

We  consider  the  opcimal  sensitivity  control  problem  formulated 
in  [Zames  1981],  but  with  plants  of  the  form  P(s)  =  A(s)+e  where 

A  and  B  are  proper  rational  functions,  and  A>0.  The  block  dieigram  in 
Figure  1  shows  the  feedback  system  we  are  considering. 


Figure  1.  Feedback  system  considered. 


The  problem  is  to  minimize  the  31  norm  of  the  weighted  closed  loop 
sensitivity,  X(s)  =  W(s)[l+P(s)F(s)]"*  over  all  stabilizing  feedbacks 
F(s). 

In  this  work  we  consider  only  weighting  functions  which  are  ratio- 

2 

nal  and  bounded  away  from  zero  at  <>.  The  reason  for  the  rationality 
assumption  is  simple  convenience,  and  the  fact  that  nothing  seems  to  be 


le  importance  of  this  was  suggested  by  Gunter  Stein. 
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lost  from  this  restriction.  If.  however,  strictly  proper  weighting 
functions  are  considered,  since  the  resulting  optimal  weighted  sensiti¬ 
vities  will  be  all-pass,  the  corresponding  unweighted  sensitivities  will 
be  unbounded  at  This  would  be  an  unacceptable  design,  were  it  really 
to  be  used.  As  above  for  the  question  of  ill-posedness,  it  could  be 
argued  that  real  implementations  will  use  proper  approximations  to  the 
optimal  compensators.  Then  the  unweighted  sensitivity  will  not  only  be 
bounded  at  Infinity,  it  will  have  value  1  there.  However,  since  the 
approximation  techniques  presented  in  the  literature,  [Zames  and  Francis 
1983]  and  [Vidyasagar  1985],  achieve  sensitivity  with  norm  arbitrarily 
close  to  the  optimal  by  approximating  the  optimal  weighted  sensitivity 
to  higher  and  higher  frequencies,  these  techniques  cause  the  unweighted 
sensitivity  to  grow  as  the  approximation  improves.  A  trade-off  would 
then  be  required.  Since  bounded  sensitivity  is  a  real  design  con¬ 
straint,  which  proper  choice  of  non-strtctly  proper  weighting  function 
accommodates,  the  use  of  such  a  weighting  function  seems  to  be  the  natu¬ 
ral  way  of  providing  the  trade-off. 

Our  solution  will  start  with  a  more  narrow  problem,  and  proceed 
stepwise  to  generalize  the  solution.  The  sequence  of  problems  to  be 
considered  is: 

A.  A=0  and  B  is  minimum  phase  and  stable,  with  sensitivity  weighting 
function  W  having  one  pole  and  one  zero.  (A=0  is  the  case  of  an 
input  delay. ) 

B.  W  is  any  stable  rational  minimum  phase  weighting  function. 

C.  B  can  have  right  half  plane  zeros. 

D.  B  can  be  unstable. 

E.  A  is  non-zero. 
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Since  the  optimal  compensator  will  in  general  be  non-proper,  we 
also  consider  the  problem 

F.  Approximation  of  the  optimum  compensator  with  proper  ones. 

As  part  of  our  solution,  we  will  examine 

G.  The  stability  properties  of  the  feedback  compensators  we  obtain, 
and 

H.  The  robustness  of  stability  of  the  resulting  closed  loop  system 
with  respect  to  perturbation  of  the  delay  interval  and  pole/zero 
locations  of  the  plant. 

In  this  report  we  discuss  problem  A.  and  the  issues  in  F  and  G  as 
applied  to  A.  Subsequent  reports  will  cover  the  other  problems. 


Derivation  of  Optimal  Sensitivity 


Formulation.  The  calculation  of  the  optimum  sensitivity  for  this 
case  was  first  presented  in  Appendix  B  of  [Flamm  1985a].  The  following 
re-derivation  is  presented  in  a  form  that  lends  itself  to  generalization 
to  the  succeeding  problems. 

The  transformation  of  the  problem  to  the  form 

taf„  IIW(s)  -  ><«(s)ll„  (1) 

Hex 

is  standard.  See.  for  example.  [Zames  and  Frtmeis  1983].  Here  W(s)  is 

the  sensitivity  weighting  function  and  is  the  inner  factor  of  the 

plant.  (We  have  assumed  that  the  plant  is  stable,  and  that  the  weight- 

ing  function  is  stable  and  minimum  phase.)  In  problem  A.  W(s)  =  and 
*sA 

'p  s  e  Note  that  we  have  normalized  the  weighting  function  so  that 
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its  zero  is  at  the  point  -1. 


Our  approach  is  based  on  [Sarason  1967],  to  which  the  reader  is 
referred  for  the  nathematical  Justification  of  our  work.  Since  we  shall 
deal  with  operators  which  are  not  compact,  we  cite  Theorem  1  on  p3Lge  179 
in  that  paper  for  the  general  case.  However,  since  we  subsequently 
establish  that  the  particular  operators  in  which  we  are  most  interested 
have  maximal  vectors,  it  is  really  Proposition  5.1  on  p.  188  that  we 
require.  The  work  in  this  report  is  further  substantially  glided  by  the 
whole  of  §7.  pp.  191-194. 

In  order  to  follow  [Sarason  1967],  we  view  ^W{s)->^<(s)H(s) j  as  an 
operator  on  The  compression  of  this  operator  to  K  =  3lf^  0  is 

equal  to  the  compression  of  W(s)  on  the  same  subspace.  Call  this  latter 
operator  T  =  The  infimum  in  (1)  cannot  be  less  than  the  operator 

norm  of  T.  Theorem  1  in  [Sarason  1967]  says  that  the  desired  infimum  is 
in  fact  equal  to  IITII. 


Solution  Technique.  In  order  to  calculate  this  norm,  we  first 

2  i.  3  '' 

characterize  K  =  (pX  )  .  Suppose  f,  h  €  L  (0,®)  smd  •  denotes  Laplace 

A  A 

transform,  g  =  Sf(g).  If  f(s)€K,  then 

0  =  <?{s),e~®^h(s)>jj2.  for  all  hear* 

=<f(l(j),e  ^‘^^h(lu)>^2 
=  <f(t),h(t-A)>j^2 

by  Parseval’s  theorem.  Therefore  f(t)  =  0  a.e.  [A,®).  Conversely,  if 
f(t)€L*  and  supp(f )C[0,A],  then  h(PX^)'‘'.  Thus 

K  =  ('<<**)■*■  =  {f:f€L’[0.®)  and  supp{f )C[0, A]> 
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where  "**"  denotes  convolution.  Accordingly,  for  f€K, 


Tf  =  i£^[u(t)-u(t-A)]J  w(t-T)f  (rjdrj ,  where  m  = 


(W). 


1-e 


-sA 


It  is  easy  to  see  =  g 


-»W(s)  L.  So 


IIT^II  = 


sup 

hesr 

-sA 


.1-e 


-sA 


rl-e 


-sA 


*W(s)(~ - *^)"2 


=  sup  ll[uj(t)w(u(t)-u(t-A))h(t)](u(t)-u(t-A))ll2 


sup  ll(io(  t)w^{  t))(u(t)-u(t-A))ll2 
supp(«)C[O.A] 

Il«ll2  =  l 


=  sup  (r|u>(t)»»«(t)|*dt)^  =  imi 

Jn 


where  V  is  the  operator  on  L*(0,A)  defined  by 


(Vf)(t)  =  U(t-T)f(T)dT. 

Jq 


(2 


A  way  to  find  this  supreimun  is  to  use  the  facts  that  IIVII*  =  I1V**VII 
(via  the  definition  of  adjoint)  and  that  p(V**V)  =  IIV**VII  since  V**V  is 
normal  [Rudin  1973.  p.  282].  Therefore  IIVII  =  pCV^V)^.  If  V  is  compact 
we  need  only  find  the  largest  eigenvalue  of  V  V.  but  this  will  not  be 
the  general  case.  However,  since  in  our  case  V  is  the  identity  plus  a 
compact  operator,  we  have  only  slightly  more  complication:  Since  we 
shall  have  V  V-I  is  compact,  we  know  from  Weyl’s  theorem  that  the  spec¬ 
trum  of  V**V  and  I  differ  only  by  eigenvalues.  [Halmos  1967.  pp.  92  & 
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295].  Therefore.  <7(V  V)  S  a(V  V-I)  U  {1}. 

Thus  we  can  find  the  norm  of  T  by  finding  the  largest  eigenvalue  of 
V*V.  We  do  this  is  two  steps: 

•  Compute  V**V, 

•  Solve  the  eigenvalue/function  problem  for  this  operator  by  trans¬ 
forming  the  eigenvalue-defining  equation  to  a  differential  equa¬ 
tion. 

We  note  that  the  operator  V  is  equivalent  to  T  via  the  Laplace 
transformation.  In  particular,  there  is  a  one-to-one  correspondence 
between  eigenvectors,  and  the  eigenvalues  are  the  same.  Furthermore, 
compactness  of  V  is  equivalent  to  compactness  of  T. 

MM  M 

Computation  of  V  V.  V  is  defined  by  <x,V  y>  =  (Vx,y>,  so  we  just 
compute: 


<Vx.y> 


y(t)‘(m(t)«x(t))dt 


Using  the  fact  that 


[We  have  used  the  fact  that 
both  ID  and  y  have  their 
support  on  [0,»).] 


iD(t-T)  =  +  «>o(t-T)  =  5(t-T)  * 

we  then  have 


fy(t)- 


(»(t)wx(t))dt 


®o(i~‘'’)x(T)dT  ♦  5j»»x(t)]dt 


=  ry(t)*[(l-/3)*e"^'*  +  x(t)] 


dt 


November  22 .  1985 


Page  10 


=  L(|  y{T)e  '"dT)(|  e^'*x{T)dT)j  - 

Jq  'Jq  0 

I  [/  y(T)e"^'^dT]*[e^'*x{t)]*dt  +  T y(t)x(t)dt 

=  [xCt)  |y(T)wo(T-t)dTdt  -  (xCt)  ry(T)u»o(T-t)dTdt 

Jq  Jq  ''0  *^0 


fx{t)[5(t)-y(t)]dt 

Jo 


=  rx(OCry(T)wo(T-t)dT  +  (6»y){t)]dt 

Jq  Jfc 

=  <  y(T)io('r-< 

''t 


-t)dT.x{t)> 


and  we  conclude 


y(T)u)(T-t)dT. 


(3) 


Solution  of  eigenvalue  problem.  We  want  to  solve 

X*f  =  vV  =  (Vo+5)»»C(Vo+6)Mf] 

=  (VoVo  +  Vo  +  Vo  +  I)(f).  (4) 


Now  let 


and  let 


y  =  Vf 


»(t-T)f(T)dT  for  t€[0,A], 


u)(T-t)y(T)dT  for  t€[0.A]. 


November  22,  1985 


Page  1 1 


•.  's. 


Since  »(t)  s  6^  +  (1-/3) *6  ,  we  can  take 

57  *‘  =  +  F 

y  =  (1-P)*xt  +  F 

1^X2=  P*X2  -  (1-P)«y 
2  *  X2  +  y 


as  a  state  space  model  for  V  V  valid  on  [O.A]. 

Bovmdary  conditions  are  given  by  y(0)  s  f(0)  and  z(A)  =  y(A).  This 
is  equivalent  to  xi(0)  =  0  and  X2(A)  =0  (if  P  ^  1:  we  see  later  that 
this  is  no  restriction). 

More  concisely. 


2  =  il-P  1]‘X  ♦  f 

0  =  [^‘(0)1 

Lxa(A)J' 


Now  we  set  z  =  X’f  in  order  to  find  the  eigenfunctions  euid  eigen¬ 
values  of  V^V.  Then 

(X*-l)f  =  [1-/3  l]-x  (e 

replaces  (5b). 

Differentiating  this  we  get 

(X*-l)f  =  [1-p  l]-i 


=  [1-D  !]•[  °1-X 

L-(i-P)*  pj 


Differentiating  again,  we  get 
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=  p^-(z-f)  +  (P*-l)f 
=  p^*(X*-l)*f  -  f  +  P*f 
=  (P*x*  -  l)-f. 


If  =  1  then  (p^-l)f  =  0  implies  f  =  0  or  p*  =  1.  This  latter 

case  implies  P  =  ±1.  and  then  our  frequency  domain  weighting  function 
S"“l 

must  be  W(s)  =  or  W(s)  =  1.  These  are  of  no  Interest  since  we  only 
allow  weighting  functions  which  are  non-constant  and  outer.  Thus  we  can 
assume  that  X^  >!  1,  and  we  get 


e!i!iLf 


Boundary  conditions  follow  from  (5c).  (6)  tuid  (7): 


(X*-l)-f(A)  =  [1-p  l]-x(A) 

=  {l-P)xt(A) 

=  (l-p)-|^e’^^^'‘^V(T)dT 

Jq 

(X*-l)f(0)  =  [-2(1-P)*  p]-x(0) 
=  P*X2(0) 

=  p.(X’-l)-f{0). 

f{A)  =  fe-P(^-"Jf(T)dT 
X*-l 


(S) 


(9a) 
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and 


f(0)  -  Pf(0). 


(9b) 


Now  it  app>ears  that  we  could  have 

(1)  3  0.  (li)  >  0  or  (iii)  2!^^  <  o. 

X*-l  X*-l 

We  consider  each  of  these  cases  in  turn. 


Case  f i 1 .  This  implies  f=0.  Then  f(t)=fc*t  +  a.  but  f(O)=0f{O) 
implies  f(t)sm(Pt+l) .  and  we  normalize  by  assuming  asl.  Then  (i)  gives 

X“P*  =  1. 


or 


X*  =  U 


The  boundary  condition  at  t  s  A  gives  us 
+  1  a 

This  implies 


X*-lJo  X*-l  Jq 


(Pi.ll(x»-U  .  r‘(Sr.l).^<iT  -  -  r*^dT  .  I*. 

Jq  0  Jq  Jq 


Pt, 
v  dr 


Ae' 


PA 


Thus  (X*-l)  =  ^*pA+l ’  substituting  X*  =  we  get 
(A+1)*P*  -  AP  -  1  s  0.  This  implies  that 


P 


_  Af>4M»CA^i}.  ^  ^  _1_ 
2*(A+1)  •  ■ 


A+r 


By  assumption  P>0,  and  so  we  must  have  P=l,  But  then  W(s) 
Thus  case  (i)  is  excluded. 


s-«-l 

s+1 


=  1. 
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Case  f 1 1 ) .  In  this  case 

f  =  e^^  +  a*e  with  k~= 

For  later  reference  we  note  this  implies 


k^-1 


(10) 


We  assume  without  loss  of  generality  that  h>0.  Then  the  boundary  condi¬ 
tion  (9b)  gives  us 


k  -  a*k  =  P*(l  +  a),  or  a  = 


k+p' 


(11) 


Boundary  condition  (9a)  gives 


Ak  -Ah 

e  -  a*e 


.  -tie. 

X’-I  ■*0 

X*-l  '>0'>  J 

1-p  -PA 

“  x’-i  *  ‘ 


^*-1  ^ 


k+p  p-k 

^(k+P)A  ^(p-h)A 


k+p 


k+P 


■)' 


using  (11).  A  little  algebra  gives  us 


and 


so 


Ah  h-0  -Ah  1-P  i  f  hA  -hA] 

*  ■  h^*  *  5^ k+p  [*  ■  *  J 

h*(e'^  -  e’'*^)  ♦  p*(e'^^  +  e~^)  »  -  e"*^) 

X*-l 
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h  *  P*coth(?t*A)  a  -  ■* 

X*-l 


Substituting  for  from  (10)  we  get 
0  s  X  ^•coth(X*A)  +  — 

P+1 

s  h  *  P*coth(X*A)  +  (k+P)*— 

P+1 

a  fe  +  p*coth(h*A)  +  (k+^)‘(—  -  1) 

P+1 


>  0 


since  k>0,  A>0.  and  therefore  coth(XA)  >  1.  This  is  a  contradiction, 
and  thus  case  (ii)  does  not  occur. 

Case  f i i  1 1 .  In  this  case 

f(t)  =  cos(<>it-4-<f)  with  u*  a  -(P  ^  — ) .  (12) 

X*-l 

Note  that  this  definition  of  w’  implies  X*  a  1  +  ^ — .  Evaluating  the 

o*+P* 

initial  condition  in  (9b)  we  see  f(0)  a  Pf(0)  implies  -usin(<p)  - 
Pcos(<f>).  implies  cos(f)^.  So  tan(<^)  a  and  therefore 


(13) 


The  initial  condition  in  (9a)  gives 


:os(uA  +  <^)  =  f (A)  a  -~P  e  P e^^cos(uT  +  <p)dT 

X*-l 

_  1~P  ^~PA  ruslnfAciH-'p)  Pco3(A<j»«^)  _  h>sin(»l  +  Pcosf<p)-i 

“  ^a_,  *•  -PA,..2^rta>  ..2.o2 


X*-l 


‘(«*+P") 


U*+P* 


_  1-0  ^(i>sinf(jA-«-«p)  -f  Pco8(A<j*«p)  _  ^-0A_usin('p)  +  Pcosf<p)j 


X*-I 


u^+P" 
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,2 


•  - rr  [cos(<f)sin(uA't^)  -  sin(<p)cos(uA+<p} 

-  e  ^^(cos{<^)sln(^)  -  sin(<p}cos(<f))] 


1-0  ^  slnfgiA) 

X*-l  (0*-Hj’)^ 


Expanding  the  left  hand  side,  we  get 

cos(uA)cos(<f}  -  sin(A<j)sin(<^) 


( l-0)sln((jA) 


(X»-1)*(«%0=')' 

Combining  this  with  the  conditions  from  (13)  we  get 


w*cos(uA)  +  0*sin(<>)A} 

X*-l 


(14) 


If  cos(<jA)=0,  (14)  implies  that  0  =  .  and  therefore  X^  = 

L*_1  P 


X^-l 


But  we  must  also  have  from  (12)  X*  =  1  +  ^  ^  .  and  so  4  =  1  + 


u*+0^ 


P 


w*+0* 


which  implies  s  0.  We  conclude  that  cos(u)=0  only  if  =  0.  Con¬ 
versely,  if  u*  «  0.  (12)  implies  X*  =  and  so  (14)  implies  wcos((i)A)  : 
0.  Then  cos(uA)  s  0  since  we  have  assumed  u^.  So  C05(uA)  s  0  is  equi- 


2n+l 


valent  to  «  =  0.  But  cos(uA)  =  0  amounts  to  wA  =  — j^ir,  an  then  from 

(13)  tan(<f)  =  -v^.  So  X*  *  i  =  ~  a  [(2n»l~)w]  eigenvalue  when 

r  -I* 

and  we  also  have 


0  =  ‘  ***  below  that  when  0 


2A 


(2n+l)wJ 


—  >  1.  It  is  the  largest  eigenvalue. 

IT 


If  cos(uA)^.  then  we  have 


w  =  -(0 


X*-l 


)tan(uA)  =  ^~^tan(uA) , 


X*-l 


Substituting  from  (12)  for  (X*-l)  we  get 
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u 


tan(uA) . 


Since  cos(uA)  ^  0  implies  jt  P,  we  conclude 


tan(uA)  =  L^*P)- 


(15) 


For  given  P  we  can  numerically  solve  this  equation  for  u,  finding  multi¬ 
ple  solutions.  From  the  definition  of  in  (12).  these  solutions  give 
us  the  eigenvalues  of  V*V  via  X*  =  1  +  ^  P  —  =  We  wish  to  pick 

the  one  among  solutions  to  the  equation  for  u  that  gives  the  largest 
eigenvalue.  We  see  that  there  are  two  cases: 

•  If  P<1  then  we  should  pick  the  smallest  solution  for  u  (usO  has 
been  excluded  by  our  consideration  of  case  (iii)). 

•  If  P>1,  we  should  pick  the  largest  solution  for  u,  but  one  can 
see  from  Fig.  1  that  there  is  no  upper  bound  on  solutions  to  (15). 
Thus  there  is  an  Infinite  sequence  of  eigenvalues  approaching  X^=l 
from  below.  (This  means  that  V  V  is  not  compact,  as  we  already 
knew. ) 

As  pointed  out  above,  the  spectrum  of  is  the  set  of  eigenvalues 
augmented  possibly  by  {1}.  Thus  for  P>1.  the  spectral  reulius  of  V  V  is 
1.  and  therefore  IIVII=1. 

If 

For  P<1,  we  also  have  l€a(V  V).  but  this  does  not  affect  the  spec¬ 
tral  radius  since  the  largest  eigenvalue  is  greater  than  1.  It  is  easy 
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to  see  that  the  eigenvalues  of  V**V  are  contained  in  Pol¬ 
and  in  (l.\~  1  for  P<1. 

Since  the  open  loop  system  has  sensitivity  of  norm  1  for  the  case 
P  >  1,  the  optimal  sensitivity  is  attained  by  H=0  in  (1),  that  is.  zero 
feedback.  In  other  words.  IIWII^  s  1  when  ^  >  1.  and  this  is  the  infimal 
sensitivity  according  to  the  above  argument.  (We  note  that  W  is  not 
inner.  For  non-compact  operators,  minimal  dilations  are  not  necessarily 
\inique  or  inner.) 

To  summarize  the  situation,  we  have  found  two  cases  of  interest: 
either  there  is  an  eigenvalue  of  V  V  equal  to  the  norm  of  T.  or  IITII=1. 
and  we  can  find  a  sequence  of  eigenvalues  of  V  V  approaching  1.  This 
will  be  used  below  to  compute  the  minimal  dilation  of  T  to  af*  in  the 
first  case.  In  the  second  case  W  itself  is  a  minimal  dilation  of  T  to 
W*. 

For  later  reference  we  note  that,  using  the  definition  of  X*  from 
(12).  (14)  Implies 

(16) 

The  solutions  of  the  equation  tan(uA)  =  oan  characterized 

graphically  as  Indicated  in  Figure  2.  (Notice  that  we  always  have  P>0.) 
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■» 


Figure  2.  Graphical  solution  of  fl5K 


corre- 


From  this  drawing  we  can  see  that  the  case  p  - 

spends  to  the  coincidence  of  vertical  asymptotes  of  y  =  tan(Ab>)  and 

y  =  If  p  s  1^1  the  asymptote  will  be  the  one  closest  to  the 

W^-P 


origin  for  tan(A(j).  Then  the  other  eigenvalues  will  correspond  to  larg¬ 
er  values  of  u.  and  thus  will  be  smaller  when  p  <  1. 


Calculation  of  optimal  sensitivity  for  S<1.  As  part  of  this  solu¬ 
tion  for  P<1  we  obtain  the  eigenfunction  f  for  the  largest  eigenvalue  of 
V  V.  This  f  is  a  maximal  vector  for  V,  and  so  we  can  follow  the  method 
in  Sarason  to  compute  the  minimal  dilation  of  T  to 

The  following  applies  to  all  eigenfunctions  f,  but  gives  us  the 
minimal  dilation  of  T  only  for  the  case  where  f  is  a  maximal  vector  for 
V. 

Let  Uo  be  a  solution  to  (15).  When  Ug  corresponds  to  a  maximal 
vector,  according  to  the  proof  of  Proposition  5.1  in  [Sarason  1967], 

A 

Tf 

T/IITII^  will  be  interpolated  by  an  inner  function  given  by  - hr.  where  • 

IITII  f 

€D 
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denotes  Laplace  transform.  In  any  case,  the  formula  below  for  will 

f 

be  valid.  First  we  calculate 

Vf  =  +  (l-/3)e'^^‘"‘^^]cos(uoT+.p)dT 

=  cos(uot-»-<f)  +  (l-P)l  e  ^^cos((iJoT+9)dT 

Jq 


=  cos(uot-t-<^)  +  (1-/3) 


sinf<Jot) 


If. 

}  "  itif) 


=  1  + 


(<jQt)*(u(  t)-u(  t~l 


=  1  + 


)  Sf  [cos  («o  t*i>)  •  (u(  t  )-u(  t-A  )  )  ] 

_ (l-/3)SfCsin(<Jat)»(u(t)-uft-A)n 


5f[ (uoCOs(uo  t  J+0sin(uo  t ) ) * (u( t )-u( t-A) ) ] 
5£[sin(«ot)-(u(t)-u(t-A))]  =  ** 

1  _.i_e“(s"^"o)A  j_^-(s+lwo)A^ 

“  2r^  ^ 

S-lUo  S+IWq 

_  2tun  -  e~^^r-(sfiuo)e^"°^f(5-t(Jn)e'^"°^~' 


2l(s’+w|) 

*  e  ^^rs«sinf(JoA)-Hjn»cosfa)nA1 


2  2 
<il  -Uo 


5£[cos(<Jot)*(u(t)-u(t-A))]  = 


,  -sA  6,  . 

1-e  „  (s-lur 


1  1  -(s-t«o)A  ,  -(s+luo)A 

=  fci-£ - +i-2 - J 


S-iUfj 


S+lUo 


s  -  e  ^  Fs^cosCcJoA)  -  <jo»sin(a>nA, 


S*+<i>o 
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Using  (IS)  and  (19)  in  (17).  we  get 


fl-g)r~<«>n-*-  e  ^^•(s»Sln<JnA-Hi)nCOS(JoA)1 


“o[~s+e  ^^(s'COSWoA-UoSlnuoA)]  +  /3[-<i>o+e  *^(s*sin£JoA+<JoCos<i»oA)] 


=  1  +  - 

-Wo(P+s)+e 


-sA 


(l-0)r-<j„+  e  ^^(s«slniJnA-KJoCOS(JoA)1 _ 

[s(uoCOS(JoA-*-/3  sinuo^)'*^o(0cos<^o^~<^oSi^oA) ] 


(20) 


When  f  is  a  maximal  vector  for  V.  (20)  is  the  optimal  weighted 


sensitivity,  and 


Tf 


=  We  can  then  compute  the  optimal  feedback 

ntcuc 


compensator  for  the  case  0  <  1  as  in  the  next  section. 

Remark:  As  A  -»  0.  X*  -*  1.  and  as  A  •+  ®.  X*  -*  4.  This  can  be 
-  max  max  p 

seen  as  follows:  From  Figure  2.  as  A  -»  0,  uo  increases.  For  large  Uo- 
CltPiVa  -»  As  A-0,  tan(Au)~Au  over  an  increasingly  larger  interval. 

<^i-P 

Thus  the  smallest  solution  of  (15)  is  given  with  increasing  accuracy  by 


u 


Au  s  US.,  or  x  Since  X’  =  we  find  X*  x  1.  When 

^  w*+j3  1+p+P^A 


A  we  see  from  Figure  2  that  uo  -*  This  gives  X^  ~ 


1_ 

P' 


Calculation  of  Optimal  Feedback  Compensator  for  B<1. 


The  optimal  weighted  sensitivity  in  (20)  can  be  written 
uo[“S+e  ®^(s*cosuoA-uosinuoA)]  +  e  ®^(s»sinwoA+«ocosUoA)  -  Uq 
«o["S+e  *^(s*cos«oA-«osinwoA)3  +  P[-«o+e  ^^(s'slnwoA+WoCOSUoA)] 
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f-wo+e  uocosuoA+e  sinuoA)  +  woe  (coswoA-wosinwoA) 


s*(-uo'*'e  wocosuoA-t-e  psinwo^)  <^oe  (^cosuoA-wosinuoA)  -  Puo 

where  wo  is  the  smallest  solution  of  (15).  We  want  to  find  the  feedback 

F  which  results  in  this  weighted  sensitivity.  We  have  X  =  W(l-FQ)  and 
_  Q  _  W-X 

F  =  l-PQ*  so  F  =  PX- 


Taking  X  =  S.  ..  g.t  ?  .  .  (s.l)D-(..p)N 

D  PN  (s*P)e"“PoN 

N  =  s(l%ie  +  1^2®  “  Wo  and  D  =  s(?i3e  *^“Wo)  +  h.4C  *  pwo 


(s+l)[s(k3e  *^“Wo)+fe4a  ^^“pwo]  *“  {s+p)[s(hie  *^-wo)+fe2e  *^"Wo 


( 1 -P ) s inuo A (uo+s* ) 


-s*[sinuoA+«ocos«oA-<i)oe*^]+s[(<Jo“/3)sinuoA-(/3+l)ci)ocosuoA+(j5+l)uoe*^]+ 


^o[uosiiKi]oA-cosuoA->-e  ] 


(l-0)sinuoA(uo+s^} 

-s^(sinuoA-KJocos(JaA)'»-s[(uo-0)slnuoA-(/3+l}uocos(>)oA]+ 


^o((>>osinuoA-cos(i>oA)+  e®  uo[s*+(P+l)s+P] 


Poe  ^^(slnuoA-HJocosuoA) 


( l-P)siwjoA(<Jo+s^) 


(JoslnuoA-wocosuoA  s  inuoA-KjocosuoA  uosinuoA-WoCOSUoA 

- p. - s  +  /3* - 

sinuoA-hJoCosuoA  sIiwoA-HjocoswoA  sinuoA-KJoCOSuoA 


+  uoe 


sA  s*+(P+l)s+P 
s  i  nuo  A-KiJo  c  o  suo  A 


Using  (16)  we  get 


(l-P)sinuoA(«o+s*) 


Po  ^(slnwoA+«ocos<i)oA)(-s*+P*)e  +  «o[s*+(P+l)s+/3]j 


=  •3ln(JoA*Po^ 
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s*+{p+l)s+p  g^-s^ 


Uo 


s*+(/3+l)s+/3 


Taking  ft  =  •^-^slnuoA  and  fa  = 
GJo 


sinuoA-HJocosuoA 


uo 


(the  expression  for 


fa  is  simplified  below)  this  is 


^  -1 
F  =  ft*Po  • 


S*+<i)o 


a^*(p+l)s-^p  j  ^  e~®^.fa _ 


a  _3 


s  +(/3+l)s+P 


which  can  be  realized  as  shown  in  Figure  3. 


Figure  3.  Realization  of  Optimal  Compensator 

Since  the  optimal  PQ  is  independent  of  the  outer  part  of  P,  Q  and 
thus  F  will  in  general  be  Improper.  It  is  necessary  to  find  a  proper  Q 
for  the  compensator  to  be  physically  implemen table. 

Proper  ApDroxlmation  of  Optimal  Feedback  Compensator. 

The  only  procedures  in  the  literature  for  this  purpose  seem  to  be 
given  in  [Zames  and  Francis  1983.  p.  591]  and  [Vidyaseigar  1985,  p.  178]. 
However  the  procedure  in  [Zames  and  Francis  1983]  requires  the 
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evaluation  of  the  term  where  3^(8)  is  the  Blaschke  product  formed 

from  plant  zeros.  In  our  case  there  is  no  Blaschke  product  involved, 
but  rather  a  singular  inner  function.  If  we  interpret  B^(s)  to  be  this 
inner  function.  B^(<»)  is  not  defined.  In  fact  this  procedure  does  not 
work  for  our  csise. 

The  procedure  in  [Vidyasagar  1985]  does  not  work  for  our  case 
either.  The  essence  of  the  difficulty  is  the  sane  as  in  the 
Zanes-Francis  procedure  —  the  inner  factor  of  the  plant  is  not  contin¬ 
uous  at  infinity. 

An  explanation  of  why  these  techniques  do  not  work  is  given  in 
[Flamm  1985b].  The  analysis  there  does  suggest  a  means  to  accomplish 
the  approximation.  A  full  discussion  will  appear  in  Part  III  of  this 
report,  but  the  graphical  argument  presented  in  [Flamn  1985b]  to  show 
why  the  other  procedures  do  not  work  motivates  the  following  approach. 

The  essential  idea  is  to  roll-off  the  ideal  Q-parameter  with  a 
transfer  function  for  which  the  Bode  magnitude  plot  has  slope  less  than 
1.  so  as  to  limit  the  phase  deviation  due  to  the  roll-off,  until  suffi¬ 
cient  attenuation  has  been  obtained.  This  this  can  be  accomplished,  for 
example,  with  a  lead-lag  network  which  approximates  such  compensation  by 
having  average  slope  magnitude  less  than  1. 

We  note  that  any  stable  Q-parameter  results  in  a  stable  closed  loop 
system,  so  that  this  roll-off  technique  preserves  stability  Just  as  the 
procedures  for  the  Isolated  right  half  plane  zero  case  do. 


November  22,  1985 


Page  26 


The  compensator  Is  of  the  form 


F  a  fi  .  Po  •  .  - i - . 

(s+P)  (s+l)+e'*^C2(s-j3) 

Since  Pq  is  minimum  phase  by  assumption,  the  problem  amounts  to  deter¬ 
mining  the  stability  of  the  term - ^ — j - .  Equivalently,  we 

(s+l)+e‘*^Ca{s-P) 

ask  whether  1+e  has  zeros  in  the  right  half  plane.  We  answer 

this  question  by  proving  that  1+e  has  finitely  many  zeros  in  the 

closed  left  half  plane,  and  then  conclude  by  appealing  to  Picard's  theo¬ 


rem  that  1+e  has  infinitely  many  zeros  in  the  right  half  plane. 

First  we  note  that  |e  >1  for  s  in  the  left  half  plane,  and 

|e  ^^1  <  1  in  the  right  half  plane.  Now  all  zeros  must  satisfy 
-Ble(sA)  ls-0  1 

e  '  ^*|s+l  “  I ■  therefore  all  closed  left  half  plane  zeros 

Is-0 I  1 

^  So  all  closed  left  half  plane  zeros  lie  on  or 

inside  the- intersection  of  the  ellipse  ||^  =  with  the  closed  left 
half  plane.  See  Figure  4. 
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Thus  all  closed  left  half  plane  zeros  lie  in  a  compact  region,  and 
we  conclude  from  analytic! cy  that  there  are  only  finitely  many  in  chac 
region.  Picard’s  theorem  tells  us  that  there  are  infinitely  many  zeros, 
so  that  we  must  conclude  that  there  are  infinitely  many  in  the  open 
right  half  plane. 

We  can  actually  get  more  detailed  information  about  the  distribu¬ 
tion  of  the  right  half  plane  zeros  without  much  more  trouble.  As  stated 
above,  the  right  half  plane  zeros  must  satisfy  j“^|  >  | .  As 

|s|  -*  •,  -»  1.  Therefore  as  |s|  -*  •.  the  zero  set  <z^>  approaches 

the  line  le”*^!  =  | .  which  is  the  same  as  Se(s)  =  and 

>a(z^)  -»  (2n-t-l)T.  Also,  since  right  half  plane  zeros  must  satisfy 
s*0 1 

<  1  (by  comparing  distances  from  P  and  from  the  point  -1).  we 
have  <  |Ca|.  Then  since  ^ 

and  we  conclude  3Je(s)  <  I  for  these  zeros. 

A 

The  fact  that  H(s)  *  (s+l)+e  fa(s-p)  has  right  half  plane  zeros 
can  also  be  obtained  from  results  on  the  distribution  of  zeros  of  entire 
functions.  (See  [Levin  1980]  Qiapter  7,  54,  p.  323,  Example  1.)  For 
our  purposes  it  is  more  convenient  to  refer  to  the  earlier  work 
[Pontryagin  19K]. 

According  to  theorem  6  in  this  latter  paper  it  is  sufficient  to 
show  that  the  function  G(y)  >  i/*cos(yA}  sln(yA)  Cs*!/  has  zeros  which 
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are  not  purely  real.  To  see  this  all  we  need  is  to  recognize  that  G 
has  infinitely  many  zeros.  (To  show  this  we  can  use  theorem  3  in  the 
same  paper.)  SettingG(y)  s  0  we  get 


-ycos(yA)  =  sin(yA)  ♦  Cs’W* 

For  real  values  of  y,  the  left  hand  side  of  this  equation  has  the  lines 
z  s  ±y  for  an  envelope,  whereas  the  right  hand  side  oscillates  with 
deviation  1  about  the  line  z  s  Since  >1,^  for  some  value  yo 

there  are  no  more  real  zeros  for  |y|  >  yo. 
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^Theorem  6  says  that  if  we  evaluate  H(s}  on  the  imaginary  axis,  and 
split  the  resulting  H(ty)  (y€]R)  into  real  and  imaginary  parts. 

H(ly)  s  F(y)  *  iG(y)  with  F(y)  and  G(y)  taking  only  real  values, 
then  necessary  and  sufficient  conditions  for  H  to  have  all  its  zeros 
in  the  left  half  plane  are  that  (1)  F  and  G  have  only  real  zeros, 
that  (ii)  these  zeros  alternate,  and  (lit)  for  at  least  one  value  of 
y.  G‘(y)F(y)  -  F'(y)G(y]  >  0.  Thus  we  need  only  show  that  G  has  a 
complex  zero  to  establish  that  H  has  a  right  half  plane  zero. 


%rom  (15),  -  — ll^cosuo^-  Using  this  in  the  defini-tion  of  f2- 


«o 


«o-P 


l-Kio 


fa  a  »cos4>oA.  Therefore  f| 
uo-P 


1+wo 


wo+1 


l-«>tan^UoA  Uo*P 


where  we 


have  used  (15)  to  substitute  for  tan(uoA).  and  simplified. 
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